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Abstract: The negative mode of the Schwarzschild black hole is central to Euclidean 
quantum gravity around hot flat space and for the Gregory-Laflamme black string in- 
stability. We analyze the eigenvalue as a function of space-time dimension A = \{d) by 
constructing two perturbative expansions: one for large d and the other for small d — 3, 
and determining as many coefficients as we are able to compute analytically. Joining the 
two expansions we obtain an interpolating rational function accurate to better than 2% 
through the whole range of dimensions including d = 4. 
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1. Introduction and Summary 

The Schwarzschild black hole is known to possess a single ofF-shell negative mode^ discov- 



In the context of hot flat space the Euclidean black hole appears as a saddle point of 
the action with the same asymptotics, namely a periodic Euclidean time direction, and is 
interpreted to represent a non-perturbative (tunneling) decay mode through the nucleation 
of a real black hole. The resulting leading correction to the spacetime energy density comes 
from the exponential of minus the action. If this non-perturbative effect is to describe an 
instability this correction should be imaginary and whether this is the case is determined 
by the subleading correction given by the square-root of the perturbations determinant. 
Therefore Gross-Perry- Yaffe argue that the consistency of this picture requires an odd 
number of negative modes, and in particular at least one mode should exist. 

The Gregory-Laflamme string instability occurs only for wavenumbers smaller than 
a certain critical wavenumber kci, which translates into a critical mass in the presence 
of a fixed compactification length. kcL has a clear physical meaning and it is directly 
determined by the negative eigenvalue, as the latter is nothing but (— fe^^). 

The negative eigenvalue was computed numerically in 4d to be ^ A(4(i) ~ 0.76 and 
extended to arbitrary dimensions in [P, where it was shown that the behavior of kciid) 
brings insight into the discovery of a critical dimension where the GL transition turns 

^This mode does not represent a physical instability of the black hole but rather of the black string 
constructed from it. 

^The data was obtained in [H] and published in [m|. 
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second order The negative eigenvalue was also computed in Q for 4 < d < 9, in Q 
for 5d, and in Harmark-Obers coordinates were used to obtain a master equation. The 
computation was generalized in [^] to boosted strings, in [11 1 to charged strings and in 
1 12] to rotating strings. Finally related non- uniform strings in various d appeared in [13|. 

While an analytic determination of X{d) is not known it is usually useful to have some 
analytic insight. In this paper we attempt for analytic control by a standard strategy in 
physics - looking at extreme limits and perturbing around them with a small parameter. 
Here, taking d as a variable allows for two limits: high d and low d where for reasons to be 
described later at low d we are interested in d ^ S"*" (approaching 3 from above). Actually 
the large d limit of this problem was already considered in but here we are able to obtain 
more results as discussed towards the end of this introduction. 

In high dimensions gravity becomes increasingly short-ranged, and this limit was 
already considered as a simplified limit for the Feynman diagrams of quantum gravity 
1 14, [l^, [l^. The idea of small, continuously varying dimensions is also familiar from the 
technique of Dimensional Regularization in field theory |17]. Moreover, 3d gravity is known 
to have an interesting and somewhat degenerate limit (see |18|,|19] and references therein) 
where gravitational waves do not exist and space-time is fiat apart from deficit-angle point 
singularities at the location of each particle. 

Recently "the optimal gauge" master eigenvalue equation for the negative mode was 
obtained in |2C, 21]. This "optimal gauge" was generalized by proving that for perturba- 
tions of any spacetime with at most one non-homogeneous coordinate the gauge freedom 
can always be eliminated 1^2| thus avoiding the need to choose a gauge fix. The lat- 
ter general theorem was already applied to a re-derivation of gravitational waves in the 
Schwarzschild background ]|2^ . 

In this paper we utilize the new master equation to analytically study the function X(d) 
for both large and small dimensions. In section ^ we set the stage by introducing the master 
equation and showing that in both the high and low d limits it can be recast in terms of 
two different variables and analytic control of A can be achieved through the technique of 
Matched Asymptotic Expansion (MAE) between the near horizon and asymptotic regions. 
In sections ^ and |^ we analyze the equations for high and low dimensions respectively. In 
section ^ we compare the analytic results of the previous two sections to new high-precision 
numerical data (summarized in appendix^) and propose an interpolating rational function 
valid over the whole range of d. The interpolation resembles the Fade approximation only 
here it is based on two Taylor expansions (rather than one) from both ends of the range 
3 < d < oo. 



Summary of results. The high d expansion is found to be 

\ = d-l + l + --- = d-A + l + ... , (1.1) 
d d 

^More insight into this critical dimension and other aspects of the problem can be gleaned from a 
hydrodynamical analogue of the Rayleigh-Plateau instability which is responsible for the creation of 
drops in a faucet. 
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where d = d — 3. See a discussion of the rehabihty of the last term of (|1.1| ) in the second 
paragraph of section |5[ The low d expansion is 



A = Cl e + C2 • • • = Cl (d - 3) + C2 - 3)2 . . . 
Cl ~ 0.71515 

C2 ~ 0.0627 , (1.2) 

where e = d = d — 3 is small, and while the leading O (e) behavior was inferred analytically 
the constants ci, C2 were determined numerically. 

The two expansions can be combined into the following rational approximation valid 
throughout the range < d < oo 

A 7 . (1,3) 

(l-ci)d + l 

The approximation was built to match both series (apart for the last term of each one), 
and it is found to uniformly agree with the numerical data to better than 2%. In particular 
by approximating ci ~ 0.7 we find the following handy formula 

A ~ d — , (1.4) 

3d + W 

and in 4d 

A(4d) « ^ = 0.769 , (1.5) 
which is in good agreement with our available numerical value A(4d) = 0.76766. 

Discussion. In general the idea to study GR in various dimensions, high and low, is 
received with considerable resistance from both laymen and professionals. In our opinion, 
one of the important motivations to consider various d's is the view that General Relativity 
is defined for all d, and hence d should be considered as its parameter. From this point 
of view it is as natural to study GR is various d's as it is to study field theories with 
matter content that differs from the standard model. We view the success of the analysis 
presented here, including the insight into 4d results as a vindication of this "dimension as 
a parameter" approach. For completeness it should be mentioned that there are at least 
two other important reasons to study GR in various d's: theoretical reasons from String 
Theory, and "phenomenological" reasons in the context of the Large Extra Dimensions and 
braneworld scenarios. 

Comparison with In the leading order high d behavior of X{d) was sought 
analytically. Here we are able to go further, getting two subleading orders as well as finding 
the low d behavior. Another difference is that the "optimal gauge" master equation was 
not available at the time of |^], but rather the transverse-traceless gauge was used, adding 
a spurious singularity to the equation and presumably obstructing the analysis in the near 
zone. 
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Open questions. It is plausible that our interpolating formulae ([L3,1.4) can be further 
refined. For instance A ~ (d — — 3)/d is interesting. Generalization to other pertur- 
bation problems is to be expected.'^ Another point is that since the Matched Asymptotic 
Expansion analysis amounts to forgetting the asymptotic region altogether in the actual 
calculation, it is possible that there is a more elegant argument to demonstrate this fact. 

Note added (v3). While our results are unchanged we added qualifying remarks on 
the assumed decoupling of the asymptotic region from the near horizon region (for both 
high and low d). 

2. Set-up 

The metric of the d-dimensional Schwarzschild black hole in standard Schwarzschild coor- 
dinates is given by 

/(r) = l-(^)'"' , (2.1) 

where vq, the Schwarzschild radius, is the location of the horizon and dd'^_2 denotes the 
metric of the round d — 2 sphere S'^"'^. Hereafter we shall use units such that rg = 1. 



The "optimal gauge" master equation for the negative mode is ||2y, ^ 

1 



drfr''-'dr + V{r) 



-0 = —Xip 



^ 2i,-m-if 

{2{d - 2y-3 - {d - l)f 

where '0 is a master field which is a gauge invariant combination of the metric perturbation 
components hu, hfiQ (see ||2^ for the full definition). Our sign conventions were chosen such 
that a negative eigenvalue is represented by a positive A. 

Note that the eigenvalue A = A(ci) is related to the critical Gregory-Laflamme wavenum- 
ber for the d + 1 dimensional black string through Q 

X = k^GL- (2.3) 

Large d is clearly an interesting limit. Numerical data suggests that in 3d A vanishes, 
d = 3 is also seen to be special from the definition of /(r) Therefore we shall attempt 

to analyze also d — > 3^. Altogether our two limits are 

• High d, namely d ^ oo where the small parameter is l/d 

• Small d, namely d ^ 3+ where the small parameter is d — 3. 
Accordingly we find it useful to denote 

d^= d - 3 . (2.4) 



*For quasi normal modes some high d behavior is aheady known 
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The appearance of the term r'^ ^ in the definition of / and in several other places 
in the master equation ( |2.2| ) suggests to define the variable 



X:=r'^ = r'^-^ . (2.5) 
The master equation turns out to have a nice representation in the X variable 

[-dxX{x -i)dx + Vx{x)] ip = -t^ ip 

V(X) -_^M±^>L_^ , (2.0) 



{2{d + l)X -{d + 2\ 



gotten by multiplying {2/1) by jd^, and here r should be considered to be a function of 
X. 

It will also be convenient to define 

f:=r-l X:=X-\, (2.7) 

such that the horizon lies at = r = X. 

In both high and low d limits we can define two regions (or zones) on the r, X axis 
namely, the near horizon region and the asymptotic region, and moreover their overlap 
will increase in size as the limit is taken. According to the method of matched asymptotic 
expansion, one solves for A = as follows. One solves for ■0 = ■0(r; A) separately in 
each region each time with only a single boundary condition - regularity at the boundary. 
The second boundary condition comes from matching in the overlap zone. Normally when 
propagating a wavefunction from two boundaries one requires V'VV' to be continuous at the 
matching point. When we have a matching region, rather than a point, can be character- 
ized as a linear combination of two basic solutions, namely we denote the general solution 
inside the overlap region (which is usually simpler) -0 = A{X) ipi{r; X) + B{X) "02 (^'S ^) where 
A, B are two arbitrary coefficients. The matching condition becomes 



A 

B 



(A) , (2. 

asymp 



where A,B from each region need to be read from the asymptotics of V'region(''; A) away 
from the boundary. This is the condition that fixes the eigenvalue A. 

Actually in the two limits under consideration, matching will not be necessary: we shall 
start the next two sections by demonstrating that in both cases it is enough to analyze 
the near horizon region and to replace the matching boundary conditions by a second 
regularity condition in the far side of the near horizon region zone. This means that both 
wave functions are localized in the near horizon region. ^ 



^Note that there are examples where this simphfication does not occur. For example, in the Id Quantum 
Mechanics problem of finding a bound state for a shallow potential, [—d^ + e V{x)]ip ~ E{e) ip the wave 
function of the bound state extends much beyond the size of the potential. 
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horizon 



Figure 1: The regions along the radial direction at high dimension. 



3. High d 

For high d we have f <C X and hence X = I occurs much before f = 1 on the r, X axis 
(see figure Accordingly we can define two regions 

• Near horizon f <^ 1 

• Asymptotic X ^ 1 

and the two regions have an overlap which increases in size as the large d limit is taken, in 
agreement with the general requirements of the method of matched asymptotic expansion. 



Asymptotic region and boundary conditions for the near horizon region. In the 

asymptotic region / — 1 and y ~ and hence the background is simply flat d-dimensional 
space-time. In this region the master equation becomes 

-^drr'^-^dri^ = -Xi^ . (3.1) 

A can be rescaled out of this equation after defining the dimensionless coordinate f := ^/Xr, 
and the solution satisfying regularity at infinity is 

where K^i^ is the modified Bessel function of the second kind of order d/2. 

Overlap region. In the overlap region both f <C 1 and X ^ 1, the A term in the equation 

becomes negligible and hence both r and X equations (2^,2.6) become 

-^dfr''~^df^ = ^ = -dxX^dx^ , 
whose general solution is 

i^ = A 

= A + B—.. (3.3) 



By examining the near horizon behavior of the Bessel function (^^) we convinced 
ourselves that to leading order 

^ = . (3.4) 

Below we shall adopt a working assumption that this boundary condition (which implies 
that the near horizon region does not receive "communications" from the asymptotic re- 
gion) is correct to all order which we compute. We view the good agreement of the analytic 
and numerical results as a confirmation of this assumption (see the second paragraph of 
section ^) . 

The eigenvalue problem. In the near horizon region X is a good variable ^ and we 
rewrite the X master equation (|2.6|) in a form convenient for the high d limit 



[-dxX{X -l)dx + V{X)] ip = -r\X) A V , (3.5) 
where the rescaled eigenvalue is 

A := A , (3,6) 

while the potential and r'^{X) are given by 

2 f 1 + 2/d] 

V{X) = — ^ 1 (3.7) 

(2X -1 + 2{X - l)/d\ 

r2 = exp logxj . (3.8) 

Zeroth order. The zeroth order equation in the near horizon region as d ^ oo is 

Lh^p = -\i^, (3.9) 
where the horizon zeroth order differential operator is given by 

L^, := -dx X{X -l)dx- ^^^^ " ^^-^^^ 

The boundary conditions for this equation are regularity both at the horizon X = {) 
and asymptotically at X ^ cxd. More specifically, near the horizon the two independent 
solutions behave as const, \og{X) and we require that the coefficient of the \og{X) piece 
vanishes. Asymptotically the two solutions behave as const, l/X and we require that the 
constant vanishes. 

The zeroth order equation ( p.9| ) has a marginally hound state, namely a normalizable 
solution with A = given simply by 

*o = ^ . (3.11) 



"An alternative change of variables which is of certain convenience is ^ := — tanh(log(/)/2) — 1/(2X — 1). 
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Indeed this solution satisfies the boundary conditions both at the horizon and at infinity 
and it has finite norm 

(^olV'o) = I (3.12) 

under the inner product which makes Lf^ Hermitian. This inner product is nothing but the 
standard norm 

/oo 
^*{X)i,iX)dX . (3.13) 

Since ipo has no nodes we also conclude that it is the ground state, and hence A = 
at 0{d^). 

We proceed by the standard perturbation theory for an eigenvalue problem. At each 
order k we find an equation of the form 

Lh V'fc = -Afc ^po + Srck , (3.14) 

where Afc, tpk are the coefficients in the 1/d expansion of the eigenvalue and eigenfunction, 
respectively 

oo oo 

A = ^ ~Xk/d' ip = Yl V'fc/d' . (3.15) 
Srck is the source term at the /c'th order which comes from the 1/d expansion of ( ^.61) and 



involves A and from lower orders. More specifically one needs to expand in ( |2.6D the near 
horizon potential, V{X), and = X"^/^ = exp(21og(X)/d). 



Equation (3.14) is solved first by taking the inner product with ipQ to find 



7 



= / / I / \ • (3-16) 

Next '4>k is solved from the inhomogeneous differential equation ( |3.14D while imposing 
regularity at the boundaries. This task is facilitated by the existence of a simple solution 
( |3.11 ) to the homogenous equation and by the use of the Wronskian. At each order there 



is a residual freedom to add to -i/'/fc any multiple of ipQ which amounts to a d-dependent 
normalization factor for ip. While it could be fixed by requiring tp^ to be orthogonal to ■0o> 
it is not necessary to do so, and we use this freedom to simplify the expression for ij^k when 
possible. 

First order. The first order correction to the potential is 
and the first source term is 

Srci = -Vi ipo . (3.18) 



^Note that the boundary condition (B.4) was used to derive this, and that a different one could have 
required to add a boundary term. 



-8- 



Hence by ( 3.16 ) the first correction to A is 

(-00 1 -00 



Solving for we find 



1 logX 



Second order. The second order correction to the potential is 

_ 8(X-1)(X + 1) 
^2 - + (2X - 1)4 (^-2^) 

and the second order source term is 

Src2 = [-V2 - 2 log X Ai] V'o + [-^i - Ai] 0i , (3.22) 

where the 2 log X term comes from expanding . 
By O) 



A, = = -1 . (3.23) 

Worn) 



Solving (331) for 1^2 we find 



1 logX + 1 , Li2(l-X) 

- (2X - 1)3 - (2X - 1)2 + 2X-1 ' ^^-^^^ 

where Li2{x) = ^^^=1 is the second poly-logarithm function. 

Third order. At this order we can compute analytically the correction to the eigenvalue 

A3 = 2 , (3.25) 

where we used 

16(X-1)2(2X + 1) 

= — — ■ ^ ^ ^ 

Here we stop the perturbative expansion as we did not obtain analytically the fourth 
order correction to the eigenvalue. 

Summarizing the high d section, we found that assuming the boundary condition ( ^^ ) 
the eigenvalue is given by 

\ = d-l+'i + --- = d-A+- + ... . (3.27) 
d d 
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asymptotic 



horizon 

Figure 2: The regions along the radial direction at low dimension. 



4. Low d 

In the low d limit the analysis runs in parallel with the previous high d analysis, but some 
features differ. We shall see that the matching problem still reduces to the near horizon 
region, but here we can determine analytically only the order of leading behavior and we 
need to resort to numerical analysis even in order to determine the corresponding leading 
coefficient. 

Since we are taking the limit d — 3 — > it will be convenient to use another notation 
for d 

e:=d = d-3. (4.1) 

Unlike the large d limit, here X <^ r, and hence r = 1 occurs much before X = 1 on 
the r, X axis (see figure ^ . Accordingly we can define two regions 

• Near horizon X <^ 1 

• Asymptotic f ^ 1 

and again the two regions have an overlap which increases in size as the low d limit is 
taken, allowing to employ the method of matched asymptotic expansion. 



Asymptotic region and boundary conditions for the near horizon region. In the 

asymptotic region X is a good coordinate, although r ^ X would continue to be useful. 
In this region 

_KM = ^!P±iL^<^«, (4.2) 



and hence the potential can be neglected (in the last inequality it should be borne in mind 
that A = C(e) as we shall see later in this section). 

On the other hand the function /(r) = 1 — r^^ = X /{I + X) varies in this region: for 
large X we have / — 1 while for small X we have / ~ X ~ e log(r). Yet, throughout the 
asymptotic region the leading expression for the solutions is 



ip = A exp(\/Ar) + B exp(— \/Ar) 



(4.3) 



Therefore regularity at infinity requires choosing A = 0. Again we shah make a working 
assumption that this boundary condition is correct to all orders which we compute and we 
view the good agreement with the numerical results as a confirmation. 



The eigenvalue problem. In the near horizon region r is a good coordinate. The master 
equation reads (^) 



1 



drfir)r^+'dr + V{r) 



^p = —A 



f 
V 



2e^ (2 + e) 
r2 (2(1 + e) (r^ - 1) + e) 



(4.4) 



For fixed r we have r"^ — 1 = exp(e log(r)) — 1 = e log(r) + . . . . Hence for small e we 
may approximate / by 

/ = e log(r) + . . . . (4.5) 

In order to have an expression for V{r) in this limit we expand the following expression 
from its denominator 



Dnv := 2(1 + e) (r^ - 1) + e = e(21og(r) + 1) + . . . 
Altogether the zeroth order eigenvalue problem is 

ip = —Xif) 



1 4 
--drr log(r) dr - ^7-^ — t \ , ^\2 
r (2 log(r) + Xy 



where in this section the rescaled eigenvalue is defined by 

A := A/e , 



(4.6) 



(4.7) 



(4.8) 



which is different from the high d definition. The scalar product which makes this operator 
Hermitian is 



(p*{r) ip{r) rdr . 



Curiously for A = there is again an exact solution 

1 - 2 log (r) 



1 + 21og(r) 



(4.9) 



(4.10) 



It is regular at both r = 1 and r = 00 but non-normalizable. For us it is important that it 
has a single node thereby indicating the existence of a unique bound state. ^ Unlike high 
d here the marginally bound state is found not to become bound for e > - which is 
fortunate, since we are expecting only a single bound state. 



*Here normalizability is not necessary. Imagine increasing A while keeping regularity at the horizon. 
During this process '4^{r) increases and the node approaches 00. The A where the node becomes infinite is 
the sought for ground state. 
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So far we have shown the existence of a single bound state whose eigenvalue scales as 
A = ci e + . . . where ci is some positive constant. In particular 



X{3d) = . (4.11) 

We were not able to solve for ci analytically. But a numerical analysis of (^]^) yields 
ci = .71515. 

One can proceed to higher orders in perturbation theory A = ci e + C2 + . . . Clearly 
this needs to be done numerically as even the zeroth order eigenvalue and eigenfunction 
are known only numerically. Doing so we obtained C2 = 0.0627. 

In summary the low d expansion is 

A = ci (d - 3) + C2 (d - 3)2 + . . . 
ci = 0.71515 

C2 = 0.0627 . (4.12) 



5. Numerical data and interpolations 

We need high-precision numerical data for X{d) in order to confirm our perturbative results 
and in order to evaluate the success of various interpolating functions based on these 
expansions. To that end we computed X{d) to 5-6 digits of precision by fairly standard 
numerical procedures and the results are collected in a table is appendix |A[ 



We first confirm that the expansions which we derived (1.1, L2|) are consistent with 
the data. The method is straightforward: to test an order k Taylor expansion around 
ro assuming order k — 1 was confirmed already, we deduct the partial sum up to order 
k — 1 from the data, divide by the small parameter to power k and seek to identify a limit 
as r — > ro (extrapolations may be used to increase the accuracy of the limit). Using this 
method we confirm the first two coefficients in ( |1.1|J1.2| ) , thereby giving us added confidence 
in the validity of our boundary conditions. The agreement with the 2/d term in the high d 
expansion is less clear. The remainder function which estimates the coefficient of the 1/d 
term is indeed very close to 2 for 12 < d < 21 but then for higher d it increases and gets as 
high as 2.5 and 3. We interpret this increase to be due to either the diminishing numerical 
accuracy for high d caused by the multiplication by (i or to a correction to the boundary 
condition ( |3.4|) at this order. 

Next we would like to exploit our control over the function at the two boundaries in 
order to construct an interpolating function which would be a reasonable approximation 
across the whole range 3 < d < oo. We shall require the function to coincide with the 
derived Taylor expansion on both side up to a prescribed order. A popular and simple 
kind of interpolating function is provided by the class of rational functions. Note that 
while the Pade approximation also employs rational functions, it is usually fitted to a 
Taylor expansion at a single point and is therefore an extrapolation which is much less 
reliable than the current interpolation. 
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It is observed that X/d is a slowly varying function. A simple possibility for a rational 
interpolating function is 

A + ai/d 



d bo + bi/d 



(5.1) 



This equation has 3 unknowns (an overall multiplication of all constants does not change 
the function) and thus can be made to fit 2 of the coefficients at high d and the single ci 
coefficient of low d. The resulting rational function is 



A{d) = d 



(1 - ci)d + ci 
{l-ci)d+l 



(5.2) 



The approximation A{d) agrees quite well with the data for X{d), as can be seen from 
figures ^J^. 
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Figure 3: The negative eigenvalue A as a function of d the space-time dimension. The blue dots 
represent the numerical data (see appendix ^) while the solid green line represents the interpolation 
Q with ci = 0.71515. 



We also attempted interpolation with rational functions of higher order, and we found 
that not only does their approximation error not improve relative to (|5.2| ) but it is actually 
even larger. 
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Figure 4: The error of interpolation in percent as a function of d (the data is represented by the 
heavy dots and the sohd curve is intended merely to guide the eye). It is seen that the maximal 
error is less than 1.7% and for d > 10 it is even less that 1%. 



A. High precision numerical values for X{d) 



In this section we present in table |l] high precision numerical results for X{d) which is 
defined by the eigenvalue problem ( |2.2D . For more details on the method of calculation see 
the Mathematica notebook p6[| . 



d 


4 


5 








A 


0.76766 


1.61015 








d 


6 


7 


8 


9 


10 


A 


2.49879 


3.41739 


4.35618 


5.30901 


6.27189 


d 


11 


12 


13 


14 


15 


A 


7.24210 


8.21779 


9.19766 


10.1808 


11.1664 


d 


16 


17 


18 


19 


20 


A 


12.1541 


13.1434 


14.1341 


15.1259 


16.1186 


d 


22 


24 


26 


28 


30 


A 


18.1062 


20.0962 


22.0879 


24.0810 


26.0750 


d 


32 


34 


36 


38 


40 


A 


28.0700 


30.0656 


32.0618 


34.0585 


36.0557 


d 


42 


44 


46 


48 


50 


A 


38.0533 


40.0515 


42.0504 


44.0505 


46.0530 



Table 1: Numerically computed high-precision negative eigenvalues A of (2.2) (in units of ) 
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